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Abstract 

In this paper we present two algorithms for the computation of a di- 
agonal form of a matrix over non-commutative Euclidean domain over 
a field with the help of Grobner bases. This can be viewed as the pre- 
processing for the computation of Jacobson normal form and also used 
for the computation of Smith normal form in the commutative case. We 
propose a general framework for handling, among other, operator alge- 
bras with rational coefficients. We employ special "polynomial" strategy 
in Ore localizations of non-commutative G-algebras and show its merits. 
In particular, for a given matrix M we provide an algorithm to compute 
U, V and D with fraction-free entries such that U MV = D holds. The 
polynomial approach allows one to obtain more precise information, than 
the rational one e. g. about singularities of the system. 

Our implementation of polynomial strategy shows very impressive per- 
formance, compared with methods, which directly use fractions. In par- 
ticular, we experience quite moderate swell of coefficients and obtain un- 
complicated transformation matrices. This shows that this method is well 
suitable for solving nontrivial practical problems. We present an imple- 
mentation of algorithms in Singular:Plural and compare it with other 
available systems. We leave questions on the algorithmic complexity of 
this algorithm open, but we stress the practical applicability of the pro- 
posed method to a bigger class of non-commutative algebras. 
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1 Introduction 

The existence and computation of normal forms of matrices over a ring is a 
fundamental mathematical question. The proof for the existence of a normal 
form is mainly constructive and can be turned into an algorithm. However, such 
a direct algorithm is not very efficient in general. Computer algebra focuses its 
attention on this kind of problems, since they are of elementary interest but of 
high complexity. 

In that sense nearly any computer algebra system is able to compute the 
Smith normal form for a matrix over a commutative principal ideal domain (Z or 
K[x] for a field K). There are many textbooks giving a theoretical background, 
like for instance [TU1[2^] . 

We present a method, which is based on Grobner bases. In [T5], there is a 
Grobner basis-based algorithm for the computation of Smith normal form of a 
matrix with entries in if [2]. Despite the fact that this approach seems to be 
folklore, we were not able to find other references. 

In this paper we consider non-commutative skew polynomial rings. Such 
rings, among other, offer the possibility to describe time varying systems in 
Systems and Control theory [32] , [17] , [18] . Many known operator algebras can 
be realized as skew polynomial rings or solvable polynomial rings [22] , some 
of them can be realized even as much easier Ore algebras [9] [SJ. However, 
general solvable polynomial rings are hard to tackle constructively (say, in a 
computer algebra system), while the class of Ore algebras of [§] IE] is indeed 
restrictive. Based on the PBW algebras [5] also known as G-algebras [33] ITS] , 
in Section [2] we propose a new class of univariate skew polynomial rings, which 
are obtained as Ore localizations of G-algebras. This framework is powerful and 
convenient at the same time. Moreover, it is more general than the class of Ore 
algebras (with defining endomorphism a being an automorphism) and allows 
algorithmic treatment of modules. In Proposition 12.21 and Theorem 12.61 several 
nice properties of such algebras (among other, these algebras are Noetherian 
domains with PBW basis) are established. We stress, that the computations 
in these algebras, especially Grobner bases for modules, are algorithmic and, 
moreover, they can be done without using explicit fractions. It is important, 
that such algebras and computations in them can be realized in any computer 
algebra system, which can handle G-algebras or polynomial Ore algebras. 

In [TT] , applications to systems of partial differential equations are shown and 
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several concrete examples are introduced. We generalize the idea, originating 
from [19], to use Grobner bases in computation of normal forms for matrices. 
The crucial improvement is introduced in Section 13.31 where we show how to 
handle the problem in a completely fraction-free polynomial framework. 

We point out advantages of the polynomial strategy and illustrate some 
of them with interesting examples in the Section where we compare our 
implementation with other available packages. In particular, we do comparisons 
with the implementation of algorithms, which use fractions directly. Notably, 
in many examples our approach delivers much more compact results with small 
coefficients. 

The non-commutative analogue to the Smith form over a principal ideal do- 
main is the Jacobson form [20] ,[10]. However, since the normal form problem is 
hard in general, we propose the notion of a weak Jacobson form, that is a diag- 
onal matrix, where the units on the diagonal will not be necessarily generated. 
Otherwise the advantage of the polynomial strategy is disturbed. Instead, we 
propose the splitting of the whole process of obtaining a (strong) normal form 
into the computation of a diagonal form and the computation of stronger di- 
agonal form from a given diagonal one. The latter, as we show in 14.41 14.51 and 
14.61 depends heavily on the domain one computes in, while the first algorithm 
is very general. 

Our implementation (of weak Jacobson and Smith forms) is realized as the 
library jacobson. lib [3D] for the computer algebra system Singular:Plural 
[TBI fTS] , which is freely available. The library has been already incorporated 
into the official distribution of Singular version 3-1-0. 

2 Algebras, Localizations and their Properties 

The framework of this paper is based on skew polynomial rings that are principal 
ideal domains. An important subclass of skew polynomial rings constitute so- 
called polynomial Ore rings. They are non-commutative rings possessing an 
endomorphism a and a cr-derivation to define the commutation rule of two 
elements, that is giving the extension from commutative polynomial ring to non- 
commutative. This kind of rings is used in analyzing the structure of analytic 
equations, like linear ordinary or partial differential equations or partial shift or 
difference equations with rational or polynomial coefficients, see Example 12.31 
The name is inspired by Oystein Ore, who introduced and studied this kind of 
rings. These rings were also studied, for instance in [SJ and [2"T] . 
Let if be a field and A be a if-algebra. Further let a : A — > A be a ring 
endomorphism. Then the map 5 : A — > A is called cr-derivation, if 5 is K- 
linear and satisfies the skew Leibniz rule 

S(ab) = a (a) 5(b) + 5(a)b for all a,b e A. 

For a cr-derivation 5 the ring A[d;a,5] consisting of all polynomials in d with 
coefficients in A with the usual addition and a product defined by the commu- 
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tation rule 



da = <j(a)d + 8(a) for all a 6 A 



is called skew polynomial ring or an Ore extension of A with d subject to 



It is easy to see, that any non-zero element a 6 A[d; a, 8] can be written as 
a = a n d n + ■ ■ ■ + aid + ao, where n £ No and <ij £ A. We call n the degree of 
a, sometimes it is also called the order of a. 

In describing i-T-algebras via finite sets of generators G and relations R, we 
write A — K(G \ R). It means that A is a factor algebra of the free associative 
algebra, generated by G modulo the two-sided ideal, generated by R. Hence yet 
another notation is A = K(G)/(R). 

Example 2.1. 

• Defining a :— id^ and 6 := we see, that (K[x±, . . . , x n ])[d; a, 8} = 
K[xi, ...,x n ,d] and K{x\,. . .,x n )[d; a, 8] = K(xi,. . .,x n )[d}. 

• Let A = K[x] for a field K of characteristic 0, a := id^-^j and S := 



is called the first polynomial Weyl algebra. 

Proposition 2.2. [5] Let A be a division ring, a : A — » A be an endomorphism 
and i? = .A [9; <r, 8] be an Ore extension with a cr-derivation (5. 
If a is injective (respectively bijective), then 

• (PID) R is a left (resp. right) principal ideal domain. 

• (Bezout's Theorem) for any non-zero a, b £ i? there exists the right (resp. 
left) greatest common divisor g r (resp. ge) of a, 6 and there exist s,t € R, 
such that g r = sa + £6 (resp. s', t', such that <7£ = as' + 6t'). 

• (ED) R is a left (resp. right) Euclidean domain. 

Hence, when a is bijective, there are left and right Euclidean division algo- 
rithms. 

In the next example we enlist some interesting skew polynomial rings (which 
are Ore algebras indeed, see 0). These rings are of great interest in applications, 
all of them can be addresses with our implementation, see Section @] 

Example 2.3. Let A = K[x), where if is a field of characteristic 0. 

• Let a := id K ^ and 8 := Then 



fi, 8. 



Wi (K) :=K[x] [d; id K[x] , — ] = K(x, d \ dx = xd + 1) 



B X (K) := A[d; id K(x) , £\ = K(x)(d \ dx = xd + 1) 



is called the first rational Weyl algebra. 
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• The first rational difference algebra is defined by 

Si := A [A; a, 6} = K(x)(A \ Ax = xA + A + 1), 
where a(p(x)) = p(x + 1) and S(p) = cr{p) — p for all p G K(x). 
Let q ^ be a unit (a parameter) in the ground field. 

. Let a(p(x)) = p(qx) and S := *(/(*)) = /( ^g 8) . Then 

Wf{K) := a, (£),] = tf(a;)(a | &c = g • + 1) 
is called the first rational g-Weyl algebra. 

• The first rational (/-difference algebra is defined by 

Q := A[d; a, 5} = K(x)(d q ■ xd + (q - l)x), 
where a(p) = p(qx) and S(p) — p(qx) —p(x). 
Indeed, we can work within the more general algebraic framework as follows. 

Let S be a multiplicatively closed set (see [27]) in a Noetherian integral 
domain A, such that g" S. S is called an Ore set in A, if for all si G S, ai G A 
there exist S2 G S,a,2 G A, such that a\S2 = S\d2- Then one can see, that 
formally (that is, allowing fractional expressions) s^ai — a2S 2 ^ 1 holds. 

Then one defines a ring of fractions or an Ore localization of A with 
respect to S to be a ring As (often denoted as S~ 1 A) together with an injective 
homomorphism <fi : A — > As, such that 

(i) for all s G S, 4>{s) is a unit in As, 

(ii) for all / G As, f — 0(s) _1 0( a ) for some a G A, s G S. 

The Ore property of S in A guarantees, that any left-sided fraction can be 
written (non- uniquely!) as a right-sided fraction. Moreover, given a\, . . . , a m G 
A and Si, . . . , s m G S, there exist a[, . . . , a' m G A and s' G S 1 , such that ajs' = 
s^a^ holds for each i. Thus there exist common right and common left multiples. 

Remark 2.4. Why such localizations are important? Among many motivating 
connections let us state the following. Given an A-module homomorphism tp : 
M — > N, where M, N are finitely generated. Then, if S~ X A exists, one has an 
induced homomorphism of 5 _1 A-modulcs S~ 1 (p : S~ 1 M — > S _1 iV. However, 
if one finds an appropriate multiplicatively closed Ore set S in A and proves 
that S -1 ^ : S~ 1 M — > S~ 1 N is not an isomorphism, it implies that M N as 
A-modules. This gives an important tool to check the isomorphy of modules. 
In contrast with common localizations of commutative ring at complements of 
prime ideals, we do not know a priori for which S we are looking for and how 
many different S should we examine. 

Note, that the question, whether two modules are isomorphic, is one of the 
fundamental questions in algebra. It is known to be not algorithmic in general, 
hence any partial algorithmic answer to this question is of big importance. 
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Definition 2.5. Let A be a quotient of the free associative algebra K(x\ x n ) 
by the two-sided ideal /, generated by the finite set {xjXi — CijXiXj — d^} for all 
1 < * < j < n, where cy € K* and dij are polynomials in xi, . . . ,x n . Without 
loss of generality [23] we can assume that dij are given in terms of standard 
monomials x^ 1 . . . A is called a G -algebra [2HH3j, if 

• for all 1 < i < j < k < n the expression CikCjk ■ dijXk — x^dij + Cjk ■ Xjdik — 
dj ■ dikXj + djkXi — CijCik ■ Xidjk reduces to zero modulo / and 

• there exists a monomial ordering -< on K[x\, . . . , x n ], such that for each i < j, 
such that dij ^ 0, lm(dy-) -< XiXj . Here, lm stands for the classical notion of 
leading monomial of a polynomial from K[x\, . . . , x n ]. 

We call an ordering on a G-algebra admissible, if it satisfies second condi- 
tion of the definition. A G-algebra A is Noetherian integral domain [24] . hence 
there exists its total two-sided ring of fractions Quot(^4) = AA\r }, which is a 
division ring (skew field). Assume that A is generated by x±, . . . , x n+ i and sup- 
pose that the set A n (A) = {A = . . . ,i n } \h< ... < i n , K(x h ,. . .,x in | /a) 
is a G-algebra} is not empty, where I\ = {xjXi — c^XiXj — dy | i, j € A, £ < j}. 

For any A = {ii, . . . , i n } 6 A„, let us define B\ to be a G-algebra, generated 
by {x h) . . .,x in }. 

Theorem 2.6. Let A be a G-algebra in variables xi, . . . , x n , d and assume that 
A = {ari, . . . , x n } S A„. Moreover, let B := B\ and B* = B \ {0}. Suppose, 
there exists an admissible monomial ordering -< on A, satisfying Xk -< d for all 
1 < k < n. Then the following holds 

• B* is multiplicatively closed Ore set in A. 

• (B*)~ 1 A (Ore localization of A with respect to B*) can be presented as 
an Ore extension of Quot(-B) by the variable d. 

Proof. Since B is a G-algebra itself, it is an integral domain, hence B* is mul- 
tiplicatively closed and does not contain zero. Since A and B are G-algebras 
and -< is an admissible ordering, for a relation dxj — CjXjd + dj with Cj G K* 
and a polynomial dj 6 A holds dj = or lm(dj) -< Xjd. Since Xj -< 9, then 
~< d 2 , hence dj is at most linear in d. Writing dj = dj -d + bj for a,j, bj G B, 
we define dj = CjXj + dj and thus we obtain a relation dxj = c'jd + bj, where 
Xj,c'j,bj e B. 

Then, by defining cr(xj) — CjXj +aj and S(xj) = bj for all 1 < j < n, we see, 
that a is an automorphism of Quot(_B). Thus an Ore extension Quot(i?)[9; a, 5] 
is indeed another presentation of (B*)~ 1 A as soon as B* is an Ore set in A. 

Since lm(dj) = lm(aj<9 + bj) -< ^9, both ^(aj) -< and lm(foj) -< x^f? 
hold. The latter implies, that there exist positive weights w and w±, . . . , w n for 
variables {d, x\, . . . ,x n }, such that for lm(aj) = x a and lm(6j) = one has 

WiCti < and J2i w iPi < Wj + w. In particular, this can be achieved by 
setting uj big enough. Then we follow the recipe from 5 and construct a block 
ordering from this setting. Consider an ordering -<g on A, which is a block 
ordering for blocks of variables {d}, {xi, . . . ,x n }. It means that d 3> Xj for all 
j, that is the variable d is bigger than any power of Xj. The second block is an 
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ordering -<b on B, for which lm(aj) -<b Xj holds. For instance, one can take 
-<b to be the restriction of -< to B. Then \m(dj) — max^ a (ajd,bj) -<g Xjd 
holds, hence -<g is admissible ordering on A. From the Proposition 28 of [12] 
(which holds for much more general situation), the existence of such a block 
ordering as -<g implies, that the set B* is an Ore set in A. □ 

Remark 2.7. Note, that by construction Ab* ■= {B*)~ 1 A is a Euclidean (prin- 
cipal ideal) domain by the Proposition 12.21 In particular, all but one variables 
are invertible (we call them also rational variables) . We call non-invertible vari- 
ables polynomial. In a more general setting, we like to present localizations of 
the type As* , where B is a sub-G-algebra of A, as a ring of solvable type [25] 
or, equivalently, as a PBW ring [5]. In the case of several polynomial variables, 
the analogue to the Theorem 12.61 seem to be much more involved. 

Example 2.8. To illustrate the Theorem 12.61 consider the difference algebra 
Si := K(x, A | Ax = xA + A + 1). Since A -< xA is a consequence of 1 -< x (we 
assume we are dealing with well-orderings only), <Si can be localized at both 
K[x}* and K[A\*. However, the algebra, associated with the operator of partial 
integration X\ := K(x, I \ Ix = xl — I 2 ) can be localized only at but not 

at -ftf[:r]*, since I 2 -< xl is a consequence of I -< x and any ordering, satisfying 
x -< I is not admissible for X\ . 

For many problems in module theory and in applications we would like to 
analyze complicated problems via localizing at big subalgebras. In the situation 
as above, we obtain non-commutative Euclidean domain as the result, hence we 
are interested in computing Jacobson form in this setting. One of the complica- 
tions, which arise in constructive handling of objects over such algebras, is quite 
hard arithmetics in the skew field. Several fundamental questions like the trans- 
formation of a left fraction into the right one (which is possible, since the Ore 
condition is satisfied), simplification of a one-sided fraction etc. require quite 
nontrivial and complex algorithms (like computation of syzygy modules and so 
on) to be used, see for instance [T|. Even in the commutative case the com- 
putations (even with one variable) over a transcendental extension by several 
generators are still nontrivial and resource-consuming for most computer alge- 
bra systems. Hence saying "ring R is a (non-commutative) Euclidean domain" 
does not automatically mean "computations in R are easy" . 

Remark 2.9. Let us come back to the justification of terminology. Usually, 
speaking on "operator algebra with polynomial coefficients" , one means that 
one works with the set of operators d\ , . . . , d m over a commutative polynomial 
ring, say, K[x\, ...,x n }. 

By saying "operator algebra with rational coefficients" one addresses an Ore 
extension of K{x\, . . . , x n ) by the operators <9j. 

It is important to mention, that K(x±, . . . , x n ) is a localization of K[x\, . . . , x n ] 
with respect to multiplicatively closed set K[x\, . . . , x n ] \ {0}. Thus it is enough 
to define an algebra with polynomial coefficients and then speak on different 
localizations of it. Therefore the notion of Ore localization reveals the origin of 
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various "rational" coefficients and allows to treat different localizations (among 
them e.g. passage to the torus K[x-l X , . . . , x^ 1 } C K(xi, . . . , x n )) uniformly. 

3 Grobner Bases in the Computation of a Diag- 
onal Form 

3.1 Yoga with Grobner Bases 

Let us give a short introduction to non-commutative Grobner basis theory, which 
has been studied by e. g. [3[22l[23]. Suppose, that there is a G-algebra over 
a field K, which is generated by x\ , . . . , x n , d, such that i?* = A* [d; a, 5] is an 
Ore extension of a G-algebra A*, generated by {xi}. By using the lower index 
*, we point out that we deal with structures, objects in which always have a 
polynomial presentation. A nice property of a G-algebra is that as a if- vector 
space it is generated by monomials of 

Mon(JZ.) = {x^-...-x^d k \ a e N'\ke N} = {x a d k \ x a £ Mon(A*),fe € N}. 

Based on a module ordering we define leading coefficient (lc), leading monomial 
(lm), leading term (It) and leading position (lpos) notions as usual. Let a := 
(0, . . . , 1, . . . , 0) be the i-th unit vector. 

In this paper we will compute Grobner basis of modules over i?„ with respect 
to an monomial module ordering POT (position-over-term), defined as follows. 
For r, s £ Mon(.R*), 

rei < sej i < j or if i = j then r < s, (1) 

and r < s with respect to an admissible well-ordering on i?„ , eliminating d, that 
is satisfying d ^> x n > ■ ■ ■ > X\ on ]?,. 

In R, a Grobner basis is computed with respect to the induced POT ordering, 
which takes only degree of d into account since Mon(i?) = {d k | k £ N}. 

We call a £ Ji, a strict left (resp. right) divisor of b £ i?* if and only 
if 3 / 6 R* such that af — b (resp. fa = b). Extending this notation to RS 
requires that both elements a,b £ i?* have the same leading position. Moreover, 
a is said to be a proper strict divisor of b, if either b = af or b = fa holds, where 
/ is not an unit in i?*. For two monomials mi, mi £ R„ we write mi < mi for 
the comparison with the fixed monomial ordering. We say that m\ divides mi, 
if each exponent of m± is not greater than the corresponding exponent of mi- 

Definition 3.1. Let M be a left submodule of R* and < be a monomial module 
ordering on i?*. A finite subset G C M is called a Grobner Basis of M with 
respect to <, if for every f £ M \ {0} there exists a g £ G, so that lm(<?) divides 
lm(/). 

A Grobner basis G is called reduced if and only if for any pair of polynomials 
h ^ / £ G, the leading monomial \m(h) does not divide any monomial of /. 
It can be shown, that a normalized (that is with leading coefficients 1) reduced 
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Grobner basis is unique for a fixed ordering. We recall the common property of 
a Grobner basis to be, in particular, a generating set. 

Remark 3.2. Let M C R p with a Grobner basis G and / S M, Define the 
submodule S of M to be generated by all s £ G such that lm(s) < lm(/). Then 

/eS. 

3.2 Working with Left and Right Modules 

Opposite algebra. In order to work with left and right modules over an as- 
sociative if-algebra A, one has to use both A and its opposite algebra A op in 
general. Recall, that A op is the same vector space as A, endowed with the op- 
posite multiplication: V a, 6 S A op , a *a°p b = b ■ a. A natural opposing map 
makes from a right (resp. left) A-module a left (resp. right) A op -module. There 
is an algorithmic procedure to set up an opposite algebra to a given G-algebra, 
see [23]. 

Involutive anti-automorphism. Alternatively, for "swapping sides" one can 
employ an anti-automorphism 9 of A , that is a i^-linear map, which obeys 
8(ab) — 9(b) 6 (a) for all a, b S A, which is involutive, that is 9 2 = id a- Often 
such an anti-automorphism is called involution. In classical operator alge- 
bras, particularly simple involutions are known [5]. Moreover, it is possible to 
determine linearly presented involution of a G-algebra via an algorithm (Levan- 
dovskyy et al., unpublished, see Singular library involut . lib [2] for an imple- 
mentation). A constructive advantage of using involution versus using opposite 
algebra lies in the fact, that one does not need to create opposite algebra and 
make to an object its opposite. Instead, we apply an involution to an object and 
remain in the same ring. One application of involution means that the object 
we deal with change its side from left to right or vice versa. 

An involution can be defined on matrices as follows. Let 9 : A — > A be an 
involution as above. We define the map 

9 : A pxq -> A qxp , M^(9(M)) T , 

where M T is the transposed matrix of M and 9(M) = \9(Mij)] for 1 < i < p 
and 1 < j < q. 

One can easily show that (9(B ■ C)) T = (9(C)) T ■ (9(B)) T for B 6 A pxq , 
C e A qxk . Applied twice, we get B ■ C back. 

Diagonalization. Let R be a if-algebra and a non-commutative Euclidean 
PID. Recall, that a matrix U G R pxp is called unimodular if and only if there 
exists U- 1 e R pxp such that UU' 1 = U^U = id pxp . Let M £ R pxq and 
assume, without loss of generality, that p > q. Then one can show, that there 
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exist unimodular matrices U € Rp x p and V € i? 9 * 9 such that 

mi 



UMV ■ 







There are several ways to prove this statement, all based on the Euclidean 
(and thus PID) property of the underlying ring. From now on, we assume that 
R is a localization of a G-algebra as in Remark 12.71 We present algorithms 
to obtain diagonal form together with unimodular transformation matrices via 
Grdbner bases. The main idea about the computation is the sequential alter- 
nation between the computation of a reduced Grdbner basis of the submodule, 
generated by, say, the rows of a matrix and acting by the involution 9 on a sub- 
module. In the PhD thesis [TH] this idea was applied to K{x\ (of course, without 
using an involution 9, which is superfluous in that case) in order to compute a 
Smith normal form. 



In the following, by rM we denote the left i?-module generated by the rows 
of a matrix M. Further on, by G{rM) we denote the reduced left Grdbner basis 
of the submodule, generated by rM with respect to the module ordering ([1]) . 

For the i-th row of a matrix M we write Mi and My stands, as usual, 
for the entry in the i-th row and j-th. column. With respect to the context 
we identify Q{rM) = {gi, ■ ■ ■ , g m } with the matrix \g\, . . . ,£f£J*. Define the 
degree of an element ^ m G R lxq to be the degree of the corresponding 
leading monomial, that is, deg(m) := deg(lm(m)). Since R is a PID, this 
degree measures the highest exponent in the variable d. Following the standard 
convention, deg(O) = — oo. Note that the elements of Q(rM) have pairwise 
distinct leading monomials, since they form a reduced Grdbner basis. In a 
reduced Grdbner basis \m{Q{uM)i) \ \m(Q(RM)j) if and only if Q{rM){ — 
G(rM)j. 

Lemma 3.3. Order a reduced Grdbner basis in such a way, that \m(Q(RM)i) < 
■■■<\m(g( R M) m ). Then 

" S(rM)i ' 
_ Q{ R M) m _ 

is a lower triangular matrix. 

Proof. Suppose the claim does not hold. Then there exists G{RM)i and Q(rM)j 
with lpos(aU-M)i) = \pos{g{ R M) ) for i < j. Thus \m(g( R M)i) = d a 'e k and 
\m(Q(RM)j) = d aj ek such that < aj. But then evidently \m(Q(RM)i) 
divides lm(^(^M)j), which is a contradiction to Q(rM) being reduced. □ 

Due to the previous lemma, we may assume without loss of generality, that 
the matrix Q{rM) is lower triangular. Since R is an integral domain, we define 
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the rank of a matrix M to be the rank of M over the field of fractions of R. 
Now, let us assume that p = q and M is of full rank, that is row and column 
ranks of M are equal to p. The non-square case will be discussed in Remark 
1571 

Lemma 3.4. Let I denote the left ideal generated by the elements in the last 
column of 6(9( R M)), that is, by 6(9( R M) pl ), . . . , 6{9{rM) pp ). Then 

1= R (9{ R e{9{RM))) pp ). 
Proof. Note, that due to Lemma [ 



G{rM) p1 



Q{rM\ 



6(G(rM) p1 ) 
6{G(rM) pp ) 



g 



G(r9(G(rM) 



S(rM) 



According to the definition of Q the left ideal generated by G(r0(G(rM) ) ) pp 
coincides with r(6(9(rM) p1 ), 6(G{rM) pp )). □ 

Now we can formulate the algorithm that yields the desired diagonal form. 

Algorithm 3.5 (Diagonalization with Grobner Bases). 
Input: M S R 9Xg of full rank, 6 involution as above. 
Output: Matrices U,V,D G R gxg , such that 



U, V are unimodular and U ■ M ■ V = Diag(ri, . . . , r g ) = D. 



M(°) <-M,U<- id 3 x 3 , V <- id gxg 
i <- 

while (M^' is not a diagonal matrix or i 1) do 
i «- i + 1 

Compute such that [/« • M^- 1 '> = 9{rM^-^) 

if (i = 2 0)^then 
V V ■ B{UM) 
else 

[/ <- C/( l ) • [/ 
end if 
end while 
return (C/,t/,M«) 

Theorem 3.6. The Algorithm 13.51 terminates and it is correct. 
That is, for M € R 9 * 9 , let denote the matrix we get after the i-th execution 
of the while loop. Then there exists an element k £ N such that Af^O [ s a 
diagonal matrix. If k is odd, then the while loop is repeated just one more 
time (define I := k + (k mod 2) in this case). The matrices U,V obtained in 
the last loop are unimodular and satisfy UMV = Diag(mi, . . . , m g ). 
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Proof. We prove the claim by induction on g, the size of the square matrix 
M. For g — 1 there is nothing to show. Using Lemma 13.41 the equality 
R {9((M^) gg )) = h((M('')i 9 , . . . , (M«) OT ) holds. Hence we get 

fl((M«) M > C R (9((M^) gg )) for alii 

Note that 9 preserves the degree. Then the previous inclusion implies by degree 
arguments that R {(M^) gg ) = R((M < ^ r+1 ' > ) gg ) for some r. Using Lemma l3~4l and 
(M^) gg y£ (since M is of full rank), we obtain that {M^) gg is a strict left 
divisor of (M^)i g for each 1 < i < g — 1. Then the definition of Q yields that 

( ^ 

M (r+i) = W : 



V ... (M^)) gg j 

, or, in a different notation, Af( r+1 ) = W © (Af( r+1 )) 5g , that is Af( r+1 ) is a 
block matrix. 

The (g — 1) x (g — 1) matrix M' can be transformed to a diagonal matrix via 
unimodular operations by induction. It remains to consider the transformation 
matrices U and V. For each i £ N, after executing the while loop i times, we 
obtain 

f MW = • C/ ( ^ 3) • • • t/W • M • 9(Uf>) ■ 6(UW) ■ ■ ■ 6(UM), if i is even 

\ = U {l -V ■ ■ ■ ■ ■ 9{M) ■ 6{TJW) ■ 9(UW) ■ ■ ■ 9(U^), if i is odd, 

which completes the proof. □ 

Remark 3.7. In order to extend Theorem l3.6l and Algorithm [33] to non-square 
and non-full rank matrices, we need to add suitable syzygies to U respectively V 
and zero rows respectively columns to the diagonal matrix, in order to maintain 
the initial size of M. For a computational solution it is sufficient to extend 
Algorithm 13.51 in the following way. Let M l 6 R sxt where either s = p,t = q 
or s = q,t = p in the i-th while loop. Instead of computing U l , satisfying 
IP ■ M 1 - 1 = GinM 1 ' 1 ), we compute Q{ R M) for the extended matrix M := 
[id s xs M 1 " 1 ]. Such M is obviously a full row rank matrix. Defining U % :— 
[Q{ R M) T ll ...,g{ R M) T s ] T and NP := [Q( R M)J +1 , . . . , G( R M)f ] T , it is easy to 
see that WM 1 ' 1 = M l . The matrix M l consists of the rows of M 4 " 1 ) and 
additional zero rows, such that M % £ R sxt . 



3.3 Polynomial Strategy 

We are given a matrix M over a non-commutative Euclidean domain R. In this 
section, we show our main approach of this chapter. We introduce a method 
that allows to execute the Algorithm 13.51 in a completely polynomial (that is, 
fraction-free) framework. The idea comes from the commutative case and was 
appeared e. g. in [14]. 
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Let A* be a G-algebra and A — Quot(A*). Moreover, let R — A[d;a, 8], 
such that i?* = A* [d; er, 8] is a G-algebra. Evidently i?* C R, since A* C A. 
Without loss of generality, we suppose that M does not contain a zero row. 
We define the degree of an element in i?* (or Rl* 9 ) to be the weighted degree 
function with weight to any generator of A* and weight 1 to d. Thus this 
weighted degree of / £ i?* coincides with the degree of / in R. Such degree is 
clearly invariant under the multiplication of elements in A*. 

Lemma 3.8. Let M £ R gx<1 . Then there exists an algorithm to compute a 
H-unimodular matrix T £ Rl Kg such that TM £ Rj xq . 

Proof. If M £ i?* 9 * 9 , there is nothing to do. Suppose that M contains elements 
with fractions. At first, we show how to bring two fractional elements a~ l b, c~ l d 
for «,c6i,, b,d £ R* to a common left denominator, cf. [T]. For any h%, hi £ 
A*, such that h\a — h 2 c, it is easy to see that 

(fcia) _1 (frib) = aT l fq l hib = a _1 6 and (hia)- 1 {h 2 d) = (/12c)" 1 (h 2 d) = c _1 d, 

hence (hia)^ 1 = a^ 1 ^ 1 = (hie) -1 is a common left denominator. Analogously 
we can compute a common left denominator for any finite set of fractions. Let 
Ta be a common left denominator of non-zero elements from the ?-th row of M, 
then TM contains no fractions. Moreover, T is a diagonal matrix with non-zero 
polynomial entries, so it is i?-unimodular. □ 

Remark 3.9. Note that the computation of compatible factors hi for 01,02 £ 
A* can be achieved by computing syzygies, since {(hi, hi) £ A\ \ h\a\ = hia 2 } 
is exactly the module Syz(a\, —ai) C A\. The factors hi for more Oj's can be 
obtained as well. 



Define M* := TM e J?f using the notation of Lemma l3~8l Then the 
relations _r„-M* C r M and ^M, = rM hold obviously. Thus whenever we 
speak about a finitely generated submodule rM C R lxq , we denote by rM* a 
presentation of rM with generators contained in In what follows, we will 
show how to find i?-unimodular matrices U £ i?* xp and V £ i?** 9 such that 



U{TM)V 







eRl 



xq 



Since the equality U(TM)V = (UT)MV holds and UT is a i?-unimodular ma- 
trix, our initial aim follows. 

As in the previous subsection, by G(r,M*) we denote the reduced left Grobner 
basis of the submodule r,M* with respect to the module ordering <* on i?*, 
which was already defined in ([1]). Unlike the rational case, the leading monomi- 
als of elements in R.I* 9 are of the form x" 1 ■ ■ ■ x^"d^ for on, £ N. 
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Remark 3.10. Using the polynomial strategy, two improvements can be ob- 
served. On the one hand, once we have mapped the matrix we work with from 
R 9Xq to R* xq , the complicated arithmetics in the skew field of fractions is not 
used anymore. The other improvement lies in the nature of construction of 
normal forms for matrices and the corresponding transformation matrices. The 
naive approach would be to apply elementary operations inclusive division by 
invertibles on the rows and columns, that is, operations from the left and from 
the right. Indeed, there are methods using different techniques like, for instance, 
p-adic arguments to calculate the invariant factors of the Smith form over in- 
tegers [26 , but this method does not help in construction of transformation 
matrices. Surely the swap from left to right has no influence in the commu- 
tative framework. But already in the rational Weyl algebra B\ (see Example 
12.31) . ~ is an unit in B\ and = — Comparing the multiplication with 
the inverse element, that is, with x, we see that dx = xd + 1 holds. Thus a 
multiplication of any polynomial containing d with the element - in the field 
of fractions causes an immediate coefficient swell. Since a normal form of a 
matrix is given modulo unimodular operations, the previous example illustrates 
the variations of possible representations. In section [4j we present nontrivial 
examples. Especially in the Example I4.10| the polynomial strategy dams up 
the coefficient increase in a very impressive way. 

On the other hand, switching to the polynomial framework changes the 
setup. The ring i?» is not a PID anymore, what was the essential property 
for the existence of a diagonal form over R. In the sequel, we show how that 
this problem can be resolved by introducing a suitable sorting condition for the 
chosen module ordering. Referring to the argumentation of Remark 13.31 yields 
the block-diagonal form with the block above. 



o 

□ 



□ 



G{r,M«) 



(2) 



□ 



Moreover, the rows with the boxed element have the smallest leading monomial 
with respect to the chosen ordering in the corresponding block. A block denotes 
all elements of same leading position in Q (n t M* ) . In Theorem l3.15l we show that 
these elements indeed generate rM, while in Lemma 13.131 we show that these 
elements provide us with additional information. However, this result requires 
some preparations. 



14 



Lemma 3.11. Let P be R or For M G P gxq of full rank and every 

1 < * < 5, define <x; := min{deg(a) | a G pM \ {0} and lpos(a) = i}. Then for 
all 1 < i < g, there exists hi G S(pM) of degree a* with lpos(/i;) = i. 

Proof. Recall that d ^> Xj for all j. Let / € pM with lpos(/) = ? and deg(/) = 
a,. Suppose that for all g G Q{pM) with leading position i, deg(g) > ai holds. 
Since Q(pM) is a Grobner basis, there exists g G Q{pM) such that lm(g) divides 
lm(/). This happens if and only if deg(g) < deg(/) (because R* is a G-algebra 
and R is an Ore PID), which yields a contradiction. □ 

The full rank assumption in the lemma guarantees the existence of on for 
each component 1 < i < g. Note, that over i?» the cardinality of {deg(a) 
a G pA/\{0} and lpos(a) = i} is more than one in general, hence there might 
be different selection strategies. We propose to select an element according to 
min << , , see Lemma 13.131 

Corollary 3.12. Lemma [3.111 and Lemma |3~31 yield 

deg(g( R M)i) = min{deg(a) | a £ R M \ {0} and lpos(a) = i}. 

Lemma 3.13. Let be the degree of the boxed entry with leading position in 
the i-th column, that is 

ai := deg(min{& | b G G{r,M*) and lpos(6) = i} ). 

Then for all h G rM with lpos(ft-) =iwe have deg(lm(/i)) > ai. 

Proof. Now suppose the claim does not hold and there is h G rM with lpos(/i) = 
i of degree smaller than ai . Using Lemma l3.8l there exists o G A* such that ah G 
r,M*. Then deg(ah) — deg(h) and \pos(ah) = i. Due to Lemma[3jTJ deg(/) > 
ai for all / G r^M* with leading position j, hence we obtain a contradiction. □ 

Corollary 3.14. Lemma T3. 131 and Corollary 13.121 imply, that Vl<i<j 

min{deg(a) | a G 7?M\{0}Alpos(a) = i} = min{deg(a) | a G r, A/*\{0}Alpos(a) = i}. 

Theorem 3.15. Let M G R 9 * 9 be of full rank. For each 1 < i < g, let us 
define 

b; L := min{fe | b G G(r, M») and lpos(6) = i}. 

Since M is of full rank, the minimum exists for each 1 < i < g. Note that the 
set {b\, . . . , b g } corresponds to the subset of all rows with a boxed entry in the 
block triangular form[2] Moreover p(bi, . . . , b g ) = rM. 

Proof. Let / G rM\{0}. Due to Corollary 13.141 there exists 1 < k < g such 
that lpos(6fc) = lpos(/) and deg(bk) < deg(/). Thus there exists an element 
Sfc G R such that deg(/ — Skbk) < deg(bk). Since / — Skbk G rM, Corollarv l3.14l 
implies that we have lpos(/ — Skbk) < lpos(/). Iterating this reduction leads to 
the remainder zero and thus / = Sibi. □ 
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Using the notation of the previous theorem, let G*(rM) := [b\, . . . , b g ] T , 
which is by definition a lower triangular matrix. In the sequel, let M £ R 9Xg 
be of full rank. Note that then obviously G*(rM) is a square matrix. 

Proposition 3.16. Suppose M £ R gxg is a full rank matrix and there is 
J7* £ R 1 * 9 such that £/»M* — G{r,M*). Let us select the indices 

{ti, . ..,t g } C {1,...,!} such that {(E/*M*) tl , . . . , (C^Af.)tJ = G*( fl M) (3) 

in the notation of Theorem l3.15l Then U := [(J7*)ti , ■ ■ ■ , {U*)t ] T is i?-unimodular 
in R g/ a an d f/M* = G*{ R M). 

Proof. The equality f/Af* = G*(rM) follows by the definition of U. Now we 
show that U is i?-unimodular. Note that r(UM*) = R G*( R M) = R M = rM* 
holds and UM* C R 9Xg D M*. Thus there exists V £ i? 9><9 such that M* = 
V(UM*). Then = id gxg and analogously UV = id gxg since M has full row 
rank. □ 

Lemma 3.17. The equality of the following left ideals holds: 

R (9(G*( R M) gl ), . . .,e(G*( R M) gg )) = R (G*(e(G*( R M)) gg ). 
Proof. Using the argumentation given in the proof of Lemma 13.41 we obtain 

R {6(G*( R M) g i), . . . , e(G*( R M) gg )) = R{g(9(G*( R M)) gg ). 
Note the module identities rG*( r M) = rQ(rM) =>• 

6{g*( r m))r = 9(g( R M)) R r g*(9(g*(rM))) - R g(e(g( R M))). 

According to the latter identity and to the fact that both g{9{g{RM)) and 
G*(9(G*(rM))) are lower triangular matrices, we obtain R{G(9(g(RM)) gg ) = 
R (G*(9(G*( R M)) gg ). □ 

Now we are ready to formulate the polynomial version of Algorithm 13.51 
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Algorithm 3.18 (Polynomial diagonalization with Grobner Bases). 

Input: M £ R gxg of full rank, 9 an involution on _R„ and 9 as above. 
Output: i?-unimodular matrices U,V,D S i?* X9 such that U ■ M ■ V = D = 
Diag(ri,...,r 9 ). 

Find T e R 9 * 9 unimodular such that TM e R 9Xg 

M(°) <r- TM, U <-T, V <— id gxg 

i<-0 

while M^> is not a diagonal matrix or i =2 1 do 
i «- i + 1 

Compute [/« so that [/« • M^" 1 ) = £?( fl . M^) e R 1 * 9 
Select t g } C {1, ...,/} as in © 

[/W <_ [(C/W) tl ,...,(f/«) t9 ] T 

M« f- 9(G*( R M)) 
if z =2 then 

V ■k-V • 9(UW) 
else 

U <r- ■ U 

end if 
end while 
return (U,V,M®) 

Remark 3.19. It is important to mention, that the matrices U, V, D (hence the 
elements ri as well) have entries from i?*, that is polynomials. However, U and 
V are only unimodular over R and, in general, they need not be unimodular 
over i?* for obvious reasons. In each presented example we will investigate the 
case, whether U or V will be unimodular over R* as well. After all, we come 
up with the Conjecture 14. 131 

Theorem 3.20. Algorithm 13.181 terminates with the correct result. 

Proof. Using Proposition 13.161 and replacing Lemma 13.41 by Lemma 13.171 in the 

□ 

Let us 



-d and 9{x) = x. Let 



proof of Theorem 13.61 provides the claim. 

Example 3.21. Suppose R = K(x)[d;\d, ■£] and R t = K[x][d;id, 
define an involution on i?* by 9(d) 

d 2 - 
<9 2 + 



M = 



a- 
<) 



1 

X 



e R z 



Evidently T = id 2X 2 and thus Af (0) := M, U = V = id 2x2 and i 
1: Since is not diagonal, go into the while loop 

+ x + 1 xd 



• i i— 1. Since 

where Q{r,M^) 
i =2 1 



-xd - 

-a 2 - 



a + x 2 - 
xd - a + x + 2 
a — x 

x 2 d 2 + 2xd 2 + 
xd 3 + d 3 + xd 2 



a + x 
a 2 + 29 + 1 
-a -1 

d 2 + 2x9 + 2d-x 2 
+ 5d 2 - xd + 39 - : 



1 



and 
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U +- 



x 2 8 2 + 2xd 2 + d 2 + 2x8 + 2d - x 2 - 1 


-xd - 8 + x 2 + x + 1 a;9 + 9 + x 



-xd 2 - 8 2 + x - 1 
1 



d — x —9 — 1 

2: Since is not diagonal, go into the while loop 



• i +- 2. 

2 =2 

r • 



Since 



1 xd 2 + d 2 - x + 1 
1 



MW = £( fi ,MW) and 



x 2 8 2 + 2xd 2 + d 2 + 2xd + 2d-x 2 -1 
1 





1 1 



(x + l)d 2 + 2d - 

3: Since i is even and is diagonal, the algorithm returns U and V. 
And indeed, the algorithm outputs the claimed result, since 
" x 2 d 2 + 2xd 2 + d 2 + 2x8 + 29 - 



UMV 



x 2 - 1 
1 



In view of Remark 13.191 let us analyze the transformation matrices for i?*- 
unimodularity. Indeed, V is such since it admits an inverse V . On the contrary, 
U is only unimodular over R and not over i?*, since U ■ Z = W and W is first 
invertible in the localization: 



V 



1 o 

-{x + l)8 2 + X- 28- 1 1 



23- 



+ 1)9 + x - 2 



2(9 -x) (3 + 1)9- 



-4x 2 - 8a; - 4 
5a; + 5 



Algorithm 13.181 can be extended to M G i? 9 x 9 along the lines already pre- 
sented in Remark 13.71 

Example 3.22. By executing the algorithm in the 1st rational shift algebra 
K(t, S | St = tS + S) on the same matrix as in the previous example, where d 
is replaced with the forward shift operator S, we obtain a diagonal form 
Diag( {t 2 + 3t + 2)S 2 + 2(t + 1)5 - t 2 - t + 2, 1) = 



-tS - S + t 2 + 2t tS + S + t + 2 
-S + t+1 S+l 



S 2 - 1 S + l 
S 2 + 1 s - t 



1 

-iS 2 - 25 2 - 25 + t 1 



As in the previous example, it turns out that V (but not U) is even R*- 
unimodular. 



4 Implementation and Examples 
4.1 Jacobson Form 

Let R be a left and right Euclidean domain. Inspired by the Smith form, we 
will focus on how to sharpen the result of the already discussed diagonal form. 
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Theorem 4.1. jTOJ Every matrix M € R 9Xq is associated to a certain 
diagonal matrix, namely Diag(mi, . . . , mj>, 0, . . . , 0) such that additionally 

Rrrii + iR C m^i? PI Rrrii (4) 

holds for alH = 1, . . . , min{g, q} — 1. 

Due to [201 Theorem 31] the elements rtii are unique up to similarity. Two 
elements and rii are called similar if and only if there exist a,b 6 i? such 
that 

arn^ = riift, i? = aR + riiR, R = Rb + Rrrii. 

Using the notation of the previous theorem, we call Diag(mi, . . . , mi, 0, . . . , 0) 
a Jacobson normal form of M. Note that Q is hard to tackle construc- 
tively in general, since it requires to work with the intersection of a left and 
a right ideal. This difficulty disappears if R has only trivial two-sided ideals, 
that is when R is simple. Then each matrix M possesses a Jacobson form 
Diag(l, . . . , 1, m M , 0, . . . , 0) with tom € R. 

Lemma 4.2. Let A* be a G-algebra, A — Quot(A*) and R — A[d;a,S\. Let 
U, V be unimodular and a, 6, c, d G R \ {0} such that 

UDmg(a,b)V = Diag(c,d). (5) 

Then deg(a) + deg(6) = deg(c) + deg(d). 

Proof. Due to ([5]) there exists a i?-module isomorphism 

(j) : R/aR © R/bR -> R/cR © ft/<iR. 

Since A is a skew field, <fi induces an ^-vector space isomorphism. Thus the 
^-dimensions of R/aR © R/bR and R/cR(B R/dR, which are nothing else that 
the sums of degrees, coincide. □ 

Of course, inductive argument implies that sums of degrees of diagonal en- 
tries of two diagonal presentation matrices of the same module are the same. 

Jacobson normal form in the 1st Weyl algebra. Let R be the rational 
Weyl algebra K(x)[d; 1, Jj], which is a simple domain. 

Lemma 4.3. Consider a,b € R with deg(a) > 0, b ^ and deg(6) > deg(a). 
Then there exists i £ {0, . . . , deg(6) — deg(a) + 1} such that a is not a strict 
right divisor of bx l . 

Proof. Suppose that for every i £ {0, . . . , deg(6) — deg(a) + 1} there exists a 
qt E R such that bx % — ^a. Let b — b n (x)d n + • • • + b\{x)d + b (x). Note, that 
for any k € N the equality d k x = xd k + kd k ~ 1 . Thus we define r% := bx — xb = 
Y^i=i bi^id 1 " 1 with deg(ri) = n — 1 < deg(6) and n ^ since deg(6) > 1. 
Since b — qoa and bx — qia, it follows that n = bx — xb = (q\ ~ qox)a, that is 
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a is a strict right divisor of r\ . By proceeding with bx 2 and so on, we obtain 
a sequence of non-zero polynomials rj, such that deg(&) > deg(ri) > . . . and a 
is a strict right divisor of r^. Since the degree of r$ decreases exactly by 1 at 
each step, after at most deg(6) — deg(a) + 1 iterations we obtain a polynomial 
of degree deg(a) — 1, which is non-zero. Such a polynomial must contain a right 
factor of degree deg(a), what is a contradiction. □ 

The Lemma (|4.3[) suggests an algorithm to compute the Jacobson form from 
a diagonal matrix over the rational Weyl algebra. Suppose M £ R 9Xq , where 
g = q = 2. The extension to g, q £ N is evident. Algorithm 13.181 returns 
unimodular matrices U, V such that UMV = Diag(mi, m-i). Without loss of 
generality, assume that deg(m,2) < deg(mi). 

1) If ?ri2 is a unit, we get the Jacobson form just by replacing U by Diag(l, ] 
Otherwise, choose an exponent i £ N (it exists by the Lemma [4.31) such 
that mix 1 — am 2 + b with dcg(6) < deg(m2) and b ^ 0. Then 



U. 



" 1 -a' 




mi 




' 1 x l ' 




mi b 


1 




TO2 




1 




m 2 



Replace U by 



1 -a 
1 



U and V by V 



1 x l 
1 



2) Now we apply Algorithm 13.181 to the matrix 
then 

Diag(m'i,m 2 ), where deg(m 2 ) < deg(m2). 



mi r 
m-2 



The result is 



Thus, by iterating 1) and 2) we compute U and V, such that UMV — Diag(l, ttim)- 

Remark 4.4. It seems to us, that the process of obtaining Jacobson normal 
form from an appropriate diagonal matrix can be generalized to any construc- 
tive simple Euclidean PID. Moreover, it can be applied even over non-simple 
domains. There, it is not guaranteed, that the result is so nice as Jacobson 
form, but the procedure above can simplify diagonal matrices. 

Example 4.5. Over the first rational shift algebra A = K(t)(s | st = ts + s) 
(which is a not a simple domain), we provide a counterexample for a statement, 
similar to 14.31 Consider the 2x2 diagonal matrix Di = Diag(s,s). Then 
the left module Mi = A 2 /A 2 Dx (it is of dimension 2 over K(t)) is annihilated 
by the two-sided ideal (s) and hence Di is not equivalent to a matrix of the 
form Z?2 = Diag(l,p). If it were so, due to the isT(£)-dimension of Mi and hence 
M 2 = A 2 /A 2 D 2 , we see that lm(p) = s 2 . Since M 2 = A 2 /A 2 Diag(l,p) = A/Ap, 
we have Ann^ M 2 = (p) . 

Since it is not equal to Ann^ Mi = (s), Mi ^ M2. Hence, unlike over the 
Weyl algebra (or a simple domain |10|). there are many different types of normal 
forms even for diagonal matrices. 
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Example 4.6. Consider the rational g-Weyl algebra, cf. 12.31 It is not simple 
since e. g. the ideal ((q — 1)9 + a; -1 } is a proper two-sided ideal. Denote the 
generator by /, then, by the same argumentation as in the previous example we 
can show, that Diag(/, /) is not equivalent to any matrix of the type Diag(l, g). 

Since little is known about normal forms of non-simple domains, this approach 
is very interesting to investigate in the future. 

Cyclic vector method. Indeed, the existence of Jacobson form in simple 
Euclidean PID is very strong result. In particular, it tells us that any finitely 
generated module is cyclic and its presentation is a principal ideal. The method 
of finding a cyclic vector in a module and obtain a left ideal, annihilating it, 
is used in D-module theory. J. Middeke in j^H] did some investigations of this 
question. 

Conjecture 4.7. We conjecture, that the Jacobson form for, say, square matrix 
M over a simple Euclidean domain R can be computed from the given diago- 
nal form in the following way. Let M — Diag(mi, . . . , m r ). Since ^deg(mi) 
is invariant of the module R rxr /M, this number can be used as a certificate 
for probabilistic approach. Namely, consider polynomials pi of degree at most 
deg(rrii) — 1 with random coefficients in A. Compute a generator c G R of the 
left annihilator ideal of a vector [pi, . . . ,p r ] T in R rxr /M. If degc = deg(mj), 
then Diag(l, . . . , 1, c) is a Jacobson form of M. Otherwise one takes another set 
of random polynomials pi and repeats the procedure. 

One needs the probabilistic estimations on the length of random coefficients like 
in [UJ. 

4.2 Examples, Applications and Comparison 

Implementations of Jacobson normal form. To the best of our knowledge, 
Jacobson normal form algorithm has been implemented in Maple by Culianez 
and Quadrat [11] , by Robertz et al. [H [8] , by Middeke [28] and by Cheng et al 

I3EIH2]. 

We could not locate the download version of the implementation of [IT] . The 
packages FFreduce [3j and Modreduce [5] are available via personal request 
to their authors. The implementation of J. Middeke [28] was, according to its 
author, merely a check of ideas and was not supposed to become a freely dis- 
tributed package for Maple. This package is able to compute in the 1st Weyl 
algebra with coefficients in a differential field. 

D. Robertz informed us, that his implementation [4] directly follows the clas- 
sical algorithm and it has not been specially optimized. Nevertheless, in what 
follows, we compare our implementation with the one in the Maple package 
Janet 4 on some nontrivial examples. This package is available to general 
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public. 



In packages by H. Cheng et al. modular (Modreduce) and fraction-free 
(FFreduce) versions of an order basis of a polynomial matrix M from an Ore 
algebra A are implemented. In particular, such a basis is used to compute the 
left nullspace of M, and indirectly the Popov form of Al. 



Examples. As we already pointed out in the introduction, behind diagonalized 
matrices and normal forms there are various application-driven motivations, see 
e. g. m. 

Example 4.8. For instance, consider a double pendulum with lengths £\ and 
£2- Thus £1,(2 and g are constants, that is non-zero elements of K (for details 
see [11], Example 3.2.2). The linearization of this problem leads to the system 



of linear partial differential equations in d 
matrix form with the matrix 



at ■ 



M 



£id 2 +g 




id 2 



which can be written in the 



Since the variable t does not appear in M, the ground ring for the diagonalization 
process can be thought as of A — Q(£i,£2 7 g)[d]. Thus, indeed one can compute 
the Smith normal form. 

Our implementation of the diagonal form of M on this example returns 



U : 



-1/9 
-1/9 





1/9 



and V ■■ 



gh -9hd 2 

gh -ghd 2 

1 hkd 2 + gh -hhd 4 - gh - 



-g^ 

ghd 2 



such that 

UMV= I . . ° 
gh- g& 2 

This result agrees with results, obtained in [11]. Note, that a purely fractional 
method (as well as coefficient normalization procedure) will return 1 instead 
of g(£\ — £2)- With our polynomial approach we obtain a polynomial matrix, 
which is useful for further investigations. In particular, in the current example 
we see, that setting £\ = £2 implies the drop of the rank of the Smith form 
from 2 to one, thus the properties of the corresponding system will change. In 
control theory one establishes quite different properties of the module in the 
non-generic case £\ = £ 2 - 

Remark 4.9. In [5S] the algorithm for finding so-called "obstructions to gener- 
icity" was derived and discussed. A lesson learned from that paper can be 
applied for an implementation of Jacobson form as follows. It is recommended 
to split the algorithm (resp. the implementation) into two parts. In the first 
part one computes a diagonal matrix, where the invertibles of the ground ring 
are not canceled artificially. The second part applies the normalization on the 
invertibles; this part is trivial to achieve. Note, that our polynomial algorithm 
allows one to keep a close track on suspicious invertibles due to this scheme. 
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Example 4.10. Over the first rational Weyl algebra Q(t)[d; id, 4f], consider 
the matrix 

9 2 8 + 1 
R = 9+1 d 3 -t 2 d 

29 + 1 9 3 + d 2 d 2 

A typical implementation of the Jacobson normal form returns the matrix 
D = Diag(g, f , 1) together with transformation matrices U : V G Q[t] [d; id, Jj] 3x3 
such that U RV — D. Below, we write down just the leading term of each matrix 
entry and moreover, we write "l.o.t." for "lower order terms" with respect to 
degree lexicographical ordering on Q[t][9; id, 4A. The implementation of the 
Algorithm l3~18l in Singular returns D = Diag(2i 2 rf 8 + 33 l.o.t., 1,1). The 
transformation matrices are 



■ 24 l.o.t. 

l 

\a 5 + 2 l.o.t. 



ltd™ 



+ 19 l.o.t. 


~\d 2 



1 + 



f 44 l.o.t. 


2 l.o.t. 



and 



V : 



2td 2 + 3 l.o.t. 
-2i9 3 + 2 l.o.t. 
td 8 + 28 l.o.t. 



2d 2 

-2d 3 +3 l.o.t. 
d 8 + 11 l.o.t. 



2<9 2 + 1 l.o.t. 
-2d 3 

a 8 + 16 l.o.t. 



In view of lXTOl V (but not U) is unimodular over = Q[t][9;id, 

Janet returns a matrix Diag(l, 1, (279936£ 14 + 14 l.o.t.)" 1 (279936i 14 9 8 
145 l.o.t.)), 



U = 



1 

(6t 2 +2 l.o.t.)- 1 (S 2 + 1 l.o.t.) (6t 2 + 2 l.o.t.)- 1 (9 3 +3 l.o.t.) (6t 2 + 2 l.o.t.)" 1 

«31 M32 U33 



where g = (559872t 14 +14 l.o.t.), u 31 = g- 1 (-279936t 14 9 9 +158 l.o.t.), 1132 = g- 1 (279936t 14 9 
182 l.o.t.), U33 = g" 1 (279936t 14 9 7 + 127 l.o.t.) . The right transformation matrix V 



I i9 6 + 15 l.o.t. (279936* 14 + 14 l.o.t. )" 1 (46656t 12 9 7 + 110 l.o.t. ) 

9 + ll.o.t. ~|9 7 + 15 l.o.t. (-1679614* 16 + 16 l.o.t. )" 1 (279936t 14 9 8 + 138 l.o.t. ) 
" 1 (6i 2 + 2 l.o.t. )- 4 (29 2 + 1 l.o.t ) 

Example 4.11. Consider the matrix from the Example 14.101 replacing d by 
S, the forward shift operator in the first rational shift algebra in t, s. Then the 
diagonal form, computed with our algorithm is Diag(i 12 S' 8 + 101 l.o.t., 1,1). No- 
tably, the leading coefficient in t factorizes completely. Transformation matrices 
are, as expected, more complicated as in the Example 14.101 U has only three 
entries of length bigger than 1; their lengths are 113, 116, 150. In the matrix V, 
the lengths of entries are 22, 11, 58, 20, 14, 60, 26, 17, 64 with degree in S up 
to 7. Coefficients, having more than 7 digits appear only in one entry, and grow 
up to 12 digits. The situation in the first rational difference algebra is similar, 
as a reader can see by computing with our implementation. We have computed 
all the examples from this paper in the shift and difference settings as well. 
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Example 4.12. Let R = Q(y, x)[d; id, ■£] and thus R* = Q[y, x][d; id, £]. The 
matrix M below comes from the system of partial differential equations. With 
our algorithm we obtain transformation matrices and a diagonal one: 

M= \y 2 d 2 + d + l 1 

xd x 2 d 2 + d + y ' 

~x 2 d 2 — d — y 1 
1 

where g — —y 2 x 2 d 4 — x 2 d 3 — x 2 d 2 — y 2 d 3 + xd + (—y 3 — l)d 2 + (—y — 1)3 — y. 

If we consider M £ ^(y, x)[d; id, Jj] 2x2 , we obtain the single example from 
[2"8] . Then the rational form of our result is exactly the result obtained in [2"8] . 
namely Diag(l, —Ss mod 2). Note, that in our method no computations with 
4x4 matrices as in [55] are needed. As demonstrated, our implementation works 
over finite fields as well. And, as before, the right transformation matrix V is 
unimodular even over i?* = Q[y, x] [8; id, -4-]. 

As we have seen, in all the examples above the right transformation matrix 
V was indeed unimodular over i?* . We observe this phenomenon for even more 
examples over Weyl and shift algebras. 

Conjecture 4.13. Let A* be a G-algebra and A = Quot(A«). Moreover, let 
R = A[d; cr, 6], such that R* = A*[d; a, 5} is a G-algebra. For a matrix M e RP xp 
there exist square matrices U, V, D with entries from i?*, such that UMV = D, 
where D is diagonal and U, V unimodular over R. If D has only one polynomial 
non-constant entry, then V can be chosen to be unimodular over R*. 

Application. Over R, the decomposition as above can be applied as follows. 
We start with a system of equations Mu> = in unknown functions uj — 
(wi, . . . , u>p). Since U and V are unimodular over R and U MV = Diag(c?n, . . . , d pp ), 
we obtain a decoupled system {duZi — 0}, where z — V~~ 1 uj, which is equivalent 
to Muj = over R. Note, that du — is possible, then one calls z$ a free 
variable of the system in the literature (e. g. in |31j). 

Let us analyze what can be done over i?*. Suppose, that V is unimodular 
over i?*. Then UMoj = ^ DV^uj = 0. However, since U is not unimodular 
over i?* , we have implication Muj = => U Mui = only. Let T be a matrix, 
such that TU = \&r, then, by a reasoning, similar to Lemma T3.8I there exists a 
diagonal matrix Q = Diag(. . . , qu, . . .) such that Q resp. QT have with entries 
from A* resp. For simplicity, assume that A* is commutative. Denote by 
S the multiplicatively closed set, generated by q, the least common multiple of 
{qu}. If S happens to be an Ore set in R*, then the localization exists and 

U will be unimodular over S~ 1 R*. Further computations happen in different 
branches: first in the generic S~ 1 R*, where by UMw = => Muj = and then 
in the case, determined by the relation q — 0. In the latter, one can apply 
the algorithm Genericity from [25], which delivers a disjoint decomposition of 
the set of zeros of q into locally closed sets Li. One can proceed with analysis 
of systems UMw = along Li and obtain special solutions on each Li. This 



M 



1 




9 





y 2 d 2 -d-l 1 







1 
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shows, that the left transformation matrix U in this setting carries essential 
information about the so-called singularities of a system. Note, that working 
over R we compute only generic information, while following the polynomial 
strategy over i?* allows us to make a complete description of the system. 

Clearly the decoupling, provided by a diagonal form, is of big importance 
for solving systems of operator equations with rational coefficients and for the 
structural analysis, performed in the algebraic system and control theory (see 
e. g. Theorem 8 of [3Tj). 

5 Conclusion and Future Work 

Indeed, this paper is a part of a general program on providing effective computa- 
tions within Ore localized G-algebras. Notably, polynomial strategy, which we 
described in details for the case of one polynomial variable, is one of the key ele- 
ments of the program. There is ongoing work on the implementation of Grdbncr 
bases for Ore localized G-algebras under the codename Singular: :LoCAPAL. 

Polynomial strategy brings us several advantages in practical computations. 
One of them is the generality of the overall approach. Namely, as soon as there 
is an implementation of Grobner bases for modules (and hence syzygies) over a 
G-algebra A, under some mild assumptions we are able to work effectively with 
Ore localization Ab* oi A with respect to a multiplicatively closed Ore set B* , 
where B is a suitable G-subalgebra of A (cf. Theorem 12. G|) . 

The question, whether direct computations with fractions of Ab' will be 
always outperformed by the polynomial strategy, is still open. Consider, for 
instance, the situation, where the input matrix M is given already with ratio- 
nal non-commutative coefficients. Then bringing M to the fraction-free form 
is already a nontrivial operation (as soon as we work with non-commutative 
algebra) , as indicated e. g. in the proof of Lemma 13.81 

In our opinion the answer to the above question depends both on the algebra 
Ab* and on the presentation matrix M. However, in general nontrivial com- 
putation directly using fractions in the algorithm might cause the appearance 
of enormous coefficients, as several examples demonstrate. We want to stress, 
that these examples have not been specially selected for this purpose; instead, 
we picked a couple of them from a bigger family of examples. In our opinion, 
this phenomenon is quite ubiquitous. 

Our implementation of the Jacobson normal form will be developed further 
to provide a user with the possibility to compute in more general algebras. At 
the moment, the stable version of the library [3U] supports first Weyl, shift and 
difference algebras. Investigation of normal forms over non-simple domains (as 
in 14.51 14.6[) is an important future task. 

Middeke [55] has reported, that the classical algorithm, computing Jacobson 
form of a matrix over the Weyl algebra over a differential field is polynomial- 
time. However, it seems to us (due to polynomial strategy approach), that the 
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subalgebra of invcrtible elements must be involved in the complexity analysis. 
Perhaps one should consider different models for studying complexity, since ex- 
perience with practical applications suggests, that the important role, played by 
the coefficient arithmetics (which is not the arithmetics over a numerical field 
anymore!) must be appropriately reflected in the overall complexity. Otherwise 
the complexity of operations over the skew field of invertible elements remains 
hidden. 

Recently, Mark Giesbrecht and George Labahn suggested the use of another 
technique from 12 lj . namely the randomization. Starting with a matrix M, one 
multiplies M with random square (hence unimodular) matrices from both sides, 
in order to reduce the number of iterations in the Algorithms 1 and 2. Some 
experiments confirm that this might be generalized to the setting of localized 
G-algebras. However, the computations become much harder in practice due 
to increased size of polynomials to deal with. This is another reason for our 
proposal to investigate the different notions of complexity of operations over 
skew fields. 
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